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ABSTRACT 

Analysis of high spatial resolution VLA images shows that the free-free emission from 
NGC 7538 IRS 1 is dominated by a collimated ionized wind. We have re-analyzed high angular resolu- 
tion VLA archive data from 6 cm to 7 mm, and measured separately the flux density from the compact 
bipolar core and the extended (1'.'5 - 3") lobes. We find that the flux density of the core is oc , 
where v is the frequency and a is 0.7. The frequency dependence of the total flux density is slightly 
steeper with a — 0.8. A massive optically thick hypercompact core with a steep density gradient 
can explain this frequency dependence, but it cannot explain the extremely broad recombination line 
velocities observed in this source. Neither can it explain why the core is bipolar rather than spherical, 
nor the observed decrease of 4% in the flux density in less than 10 years. An ionized wind modulated 
by accretion is expected to vary, because the accretion flow from the surrounding cloud will vary over 
time. BIMA and CARMA continuum observations at 3 mm show that the free-free emission still 
dominates at 3 mm. IICO+ J = 1 ^ observations combined with FCRAO single dish data show a 
clear inverse P Cygni profile towards IRS 1. These observations confirm that IRS 1 is heavily accreting 
with an accretion rate ~2 10""* MQ/yr. 

Subject headings: accretion, accretion disks - H II regions - stars: early-type - stars: formation 



1. INTRODUCTION 

NGC 7538 IRSJ__was first detected in the radio 
at 5 GHz by iMartinI (|1973l ). who found three com- 
pact radio sources at the SE edge of the large (~ 
4') H II r egion NGC 7538 , whic h is at a distance 
2.65 kpc (jMoscadelh et all Hool). The brightest 
of the three, source B, has later become known 
as IRS 1 (|Wvnn- Williams. Becklin. fc Neugebaueil 
Il974[ ). The far-infrared luminosity of the three 
sourc es is ^^1.9 10^ L^ti, completely d ominated by 
IRS l (|Hackwell. Grasdalen fc Gehrj [1981 . IRS 1 was 
first r esolved with the VLA at 14.9 GHz by ICampbelil 
(|1984'). who showed that it has a compact (~0'.'2) 
bipolar N-S core with faint extended fan-shaped lobes, 
suggesting an ionized outflow. This is an extremely 
well-studied source with numerous masers, a prominent 
molecular outflow and extremely broad hydrogen re- 
combination lines in dicating substantial mass motion 
of th e ionized gas llGaume et al.l Il995t ISewilo et al.l 
12004 iKeto et al.ll2008D. Alth ough it has been modeled 
as an ionized jet fRevnolds '1986") or a photo-ionized 
accretion disk (Lugo, Lizano & Garay 2004), the source 
is usually referred to as an Ultra Compact or Hyper 
Compact H II region with a turnover at ~ 100 GHz. 
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Since we have carried extensive observations of IRS 1 
with BIMA and CARMA, the desire to determine the 
contribution from dust emission at mm-wavelengths led 
us to investigate where the free-free emission from IRS 1 
becomes optically thin. Such a determination is not 
possible from published data, because IRS 1 has been 
observed with a variety of VLA configurations, some of 
which resolve out the extended emission and others may 
only quote total flux. This prevented us from separately 
estimating the flux for the compact bipolar core and 
the extend ed lobes. Furthermore the flux densit y varies 
with time (jFranco-Hernandez fc Rodriguedl2004D . which 
makes it difficult to derive an accurate spectral index. 

2. VLA ARCHIVE DATA 

We have retrieved and re-analyzed a few key obser- 
vations of NGC 7538 IRS 1 from the VLA data archive. 
These are all long integration observations with very high 
angular resolution, good uv-coverage, and high image fi- 
delity. Table [T] gives the observing dates, synthesized 
beam width, sensitivity, the integrated flux of the bipo- 
lar core, and the observed total flux. We included a 
BnA arra y observation obtained by us at X-band (8.5 
GHz) (San dell. Goss fc W right 2005), because it is more 
contemporary with the high frequency data sets and pro- 
vides an additional data point at a frequency where there 
are no other high angular resolution VLA data available. 
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Fig. 1.— C-band (4.8 GHz) image of the IRSl - 3 field. The 
beam size is plotted in the bottom right corner of the image. The 
first contour is at 3-cr level, the next contour 3 times higher and 
from there on we plot 7 logarithmic contours. An outline of the jet 
from IRS 1, plotted in red, shows the optically thick inner part of 
the jet and where the jet expands and becomes optically thin. Fig[2] 
zooms in on the inner optically thick part of the jet and illustrates 
how the jet becomes smaller when we go up in frequency due to 
the outer portion of the jet becoming optically thin. 



Fig [T] shows the IRS 1 - 3 field at 4.8 GHz, while Fig [2] 
shows the morphology of the compact bipolar core as a 
function of frequency. 

3. ANALYSIS OF VLA DATA 

We determined flux densities of the compact bipolar 
core with an accuracy of a few percent. The results are 
given in Table [TJ Table [1] also gives total flux densities, 
obtained by integrating over the whole area where we 
detect emission from IRS 1. At 4.9 GHz the total linear 
extent of IRS 1 is ~ 7", while it is ~ 1'.'4 at 43.4 GHz. 
Even though we can reliably determine the flux density 
of the bipolar core at 43.4 GHz, some of the faint ex- 
tended emission is filtered out; therefore the total flux is 
underestimated. At 4.9 GHz the bipolar core has a total 
length of ~ I'.'l FWHM (FuU Width Half Maximum), 
with a separation between the two peaks of 0'.'68, at 14.9 
GHz the size is 0'.'44 with a peak separation of 0'.'24 (Ta- 
blc[T]). At 43.4 GHz we can still fit the core with a double 
Gaussian with a peak separation of ^ 0'.'12, while the lin- 
ear size (length) has a FWHM of 0'.'20. Since one of our 
data sets ( X-band, 8.4 GHz) has insufficient resolution 
to resolve the two bipolar cores, we plot the linear size 
(length) as a function of wavelength in Fig. [31 although 
fitting to the lobe separation gives virtually identical re- 
sults, but with a larger uncertainty. We find the size to 
vary with frequency as ;/-0-8±o o3^ while the flux density 
increases with frequency as !/0 'i'±o o5^ ^pj^g frequency de- 
pendence of the total flux is more difhcult to estimate, 
since some VLA configurations are not very sensitive to 
faint extended emission and will therefore underestimate 
the total flux. If we chose data sets which give a reliable 
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Fig. 2.— VLA images at 4.86, 14.94, 22.37, & 43.37 GHz showing 
the bipolar core surrounding IRS 1. Beam sizes are given in Table 
[T] The C-band image (Fig. [TJ shows faint extended lobes extending 
several arcseconds to the north and south of IRS 1. These images 
show that the core shrinks in size with increasing frequency. 

flux for IRS 2^ , which is a spherical H II region with a 
size of ~ 8", we find that the total flux from IRS 1 has 
a slightly steeper spectral index, 0.8 ± 0.03 determined 
from data sets from 4.8 GHz to 49 GHz. The fit to the to- 
tal flux densities for IRS 1 is shown as a dotted line in Fig. 
[3l Such a shallow frequency dependence would require a 
very steep density gradient in the ionized gas, if the emis- 
sion originates from an H II region ionized by a central 
0-star , although it is not impossible, see e.g. (Keto et all 
(|2008f ). A high-density H II region with a steep density 
gradient should be spherical, not bipolar with a dark 
central lane as we observe in IRS 1. Neither can a steep 
density gradient model explain the extremely broad re- 
combination line s observed in IRS 1 ([Gaume eTall [19951 : 
iKeto et all l2008l ). All these characteristics can be ex- 
plained, if t he emissio n originates in an ionized stellar 
wind or jet ([Revnoldsl [19861 . The inner part of the jet 
is optically thick and further out, where the jet expands, 
the emission becomes optically thin. At higher frequen- 
cies the outer part of the jet becomes optically thin and 
therefore appears shorter and more coUimated, exactly 
what we see in IRS 1 , see Fig [2] For a uniformly expand- 
ing spherical wind, the size of the source varies as a func- 
tion of frequ ency as t^"*^'^, while th e flux density density 
goes as j/"-^ l|Panagia fc Fellil[l975h . Since we resolve the 
emission from IRS 1, we know that it is not spherical, but 
instead it appears to originate in an initially coUimated 
bipolar jet (opening angle < 30°) approximately aligned 
with the molecular outfl ow from IRS 1 , se e Section[4l For 
a coUimated ionized jet, iRevnold^ (|1986l ) showed that it 
can have a spectral index anywhere between 2 and —0.1, 
depending on gradients in jet width, velocity, degree of 
ionization, and temperature. 

4. OUTFLOW AND ACCRETION 

Although the molecular outflow from IRS 1 in early 
studies was found to be quite compact and going from 
SE to NW ([Scoville et al.lll986l ). we find that the molec- 
ular outflow is very extended (> 4') and approximately 
aligned with the free-free jet from analyzing large mo- 
saics obtained with CARMA in ^^CO and ^^CO J = 
1 ^ 0, and HCO+ J = 1 0. For HCO+ wc filled 
in missing zero spacing with fully sampled single dish 
maps obtained with FCRAO. A thorough discussion of 
the molecular outflow is beyond the scope of this paper, 
and will be discussed in a forthcoming paper (Corder et 
al. 2009, in prep). Here we therefore summarize some of 
the main results. The GARMA mosaics show that the 

^ The emission from IRS 2 is already optically thin at 4.8 GHz 
(C-band) with an integrated flux of 1440 mjy. 



NGC,7538IRS1 - an ionized jet 



3 



TABLE 1 

Observational parameters and results of VLA archive data of IRS 1 



Frequency 


Synthesized beam 


rms 


Size 




flux density 


flux density 


Observing dates 








bipolar core 


bipolar core 


total 




[GHz] 


" X " pa = ° 


[mjy beam~^] 


" X " pa 


O 


[mJy] 


[mJy] 




4.86 


0.43 X 0.37 + 5.4° 


0.30 


1.06 X 0.44 


-9.6 


94.3 ± 4.8 


122.5 


1984:1123 


8.46 


1.21 X 0.47 - 0.5° 


0.05 


0.76 X 0.30 


-10.9 


135.7 ± 13.6=^ 


148.1=^ 


2003:1014 


14.94 


0.14 X 0.12 +47.6° 


0.017 


0.44 X 0.20 


-0.3 


185.9 ± 5.6 


248.7 


2006:0226,0512'' 


14.91 


0.14 X 0.10 +13.7° 


0.081 


0.50 X 0.20 


-1.2 


165.2 ± 13.3 


291.0 


1994:1123 


22.37 


0.08 X 0.08 -43.0° 


0.017 


0.34 X 0.17 


+3.1 


286.5 ± 11.5 


370.0 


1994:1123 


43.37 


0.14 X 0.12 +14.5° 


0.17 


0.20 X 0.14 


-1.0 


429.7 ± 12.9 


473.7 


2006:0706,0914^ 



^ Insufficient spatial resolution, blends in with IRS 2. Unreliable total flux.^ Average of two data-sets at different epoehs during the 
same year. 
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Fig. 3. — Top panel: A least squares power-law fit to the length 
of the optically thick jet as a function of frequency i/ gives 6(arcsec) 
= 4.13 u~^-^ where w is in GHz. Bottom panel: Least squares fit to 
VLA data (open squares) of the compact bipolar core surrounding 
IRSl plotted as a solid line, and to total fluxes (open triangles) 
plotted with a dotted line. In the figure we also show mm data 
points obtained with BIMA and CARMA as well as interferome- 
ter data from the literature, but they are not used in any of the 
fits. This plot clearly demonstrates that the free-free emission still 
dominates the emission in high angular resolution interferometer 
data even at 1.3 mm. 

outflow is very large and has a position angle (pa) of ^ 
—20°, which is similar to the orientation of the coUimated 
inner part of the ionized jet driven by IRS 1 on smaller 
angular scales. The outflow extends several arcminutes 
to the north and starts as a wide angle limb brightened 
flow. To the south the outflow is difficult to trace, be- 
cause of several other outflows in the giant molecular core 
in which IRS 1 is embedded, but appears more collimated 
with a linear extent of ~ 2'. Analysis of Spitzer IRAC 
archive data show that the outflow may be even larger 
to the south. The 5.8 and 8.0 /im IRAC images show a 
jet-like feature projecting back towards IRS 1, extending 
to ~ 3.5' (2.6 pc) from IRS 1, which is much further out 
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Fig. 4. — Inverse P Cygni profile in HCO+ J = 1 — > o bserved 
with CARMA supplemented with FCRAO single dish data (ICorderi 
120081) . Note that the spectrum is extracted from a continuum 
subtracted spectral line cube, hence the flux density appears to 
go below zero. The angular resolution is 4'.'5. The vertical dotted 
line shows the systemic velocity of IRS 1. 

than what was covered by the CARMA mosaics (jCordeil 
l2008f ). However, the CARMA observations also confirm 
that IRS 1 is still heavily accreting. In HCO+ J = 1 ^ 
we observe a clear inverse P Cygni profile towards the 
strong continuum emission from IRS 1 (Fig 2]). The ab- 
sorption is all red-shifted and has two velocity compo- 
nents, which could mean that the accretion ac tivity is 
yaryin g with time (episodic) . From these data ICordeil 
(f200g) estimates ~ 6.8 Mq of infalling gas. With the 
assumption that the ac cretion time is similar to the free- 
fall time, ^ 30,000 vr. iCordeil finds an accretion rate ~ 
2 10~^ Mq/yt, which will block most of the uv photons 
from the central 0-star, allowing them to escape only at 
the polar regions. 



5. MILLIMETER DATA ■ 



WHERE IS THE ACCRETION 
DISK? 



In Fig. [3] we also plotted flux densities at 3 and 1 mm 
from the literature (|Lugo. Lizano &: Garavl [20041 ) . sup- 
plemented with our own results at 3 mm from BIMA 
(Wrig ht et al. 2009, in prep) and CARMA (jCordeii 
|2008| ). All the observed flux densities at 3 mm or even at 
1.3 mm can be explained by free-free emission, with at 
most a marginal excess from dust emission. High angu- 
lar resolution CARMA continuum observations at 91.4 
and 108.1 GHz conflrm what is predicted from the fit 
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to total flux densities in Fig [3l i.e. that the free-free 
emission dominates at 3 mm. These CARMA observa- 
tions resolve IRS 1 with a size of O'.'S x <0'.'3 pa — 1° and 
0'.'4 X <0'.'l pa 20° at 94 and 109 GHz, respectively, i.e. 
the 3 mm emission is aligned with the free-free emission, 
and not with an accretion disk, which is expected to be 
perpendicular to the jet. Since we detect a coUimated 
free-free jet and since there is a strong accretion flow to- 
wards IRS 1 (Section [U, IRS 1 must be surrounded by 
an accretion disk. The morphology of the free-free emis- 
sion suggests that the disk should be almost edge-on, and 
centered halfway between the northern and the southern 
peak of the optically thick inner part of the free-free jet. 
However, if such an edge-on accretion disk is thin, i.e. the 
height of the disk is small relative to its diameter, it may 
not provide much surface area, and is therefore difflcult 
to detect at 3 mm. At 3 mm IRS 1 and IRS 2 are very 
strong in the free-free, and even the free-free emission 
from IRS 3 cannot be ignored, which makes it very hard 
to detect dust emission from the accretion disk. Since 
the spectral index for dust is '--^ 3 - 3.5, the dust emission 
may start to dominate at frequencies above 300 GHz, 
observable only by SMA, but even at 1.3 mm we may 
have a much better chance to detect the accretion disk. 
There is some evidence for such a disk in 450 and 350 //m 
conti nuum images obtained at JCMT (Sandell & Sievers 
|2004[ ). Until higher spatial resolution interferometer im- 
ages are available, it is not clear whether this emission 
originates in a disk, from the surrounding in-falling enve- 
lope, or from a superposition of several nearby sources. 

6. DISCUSSION AND CONCLUSIONS 

Free-free emission from ionized jets is very common in 
low mass protostars (Class and Class I sources), which 
all appear to drive molecular outflows. This situation is 
especially true for young stars, which often have rather 
well colhmated outflows (jRodriguez. Anglada fc Curiel 
|1999[ ). Jets may a lso be common in yo ung early B-type 
high mass stars ()Gibb fc Hoard I2007D . although they 
have not been studied as well as low-mass protostars, 
because they are short lived, more distant and harder to 
identify. Although there have not been any detections 
of jets in young O stars, i.e. stars with a luminosity 
of > 10^ L0, such objects are likely to exist. Broad 
radio recombination line (RRL) objects, many of which 
are class ifled as Ult ra Compact or Hype r Compact H II 
regions (iJaffe fc Ma rti'n-Pintado 1999; Kurtz fc Francol 
I2002t ISewilo et all 12004 : .Keto et al., ■2001 ) offer a logi- 
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cal starting point, because they show evidence for sub- 
stantial mass motions, which would be readily explained 
if the recombination line emission originates in a jet. 
Several of them also s how evidence for accretion. I n 
the sample analyzed bv 'Ja ffe fc Martm-Pintadol (jl999"), 
which partly overlap s with the sources discussed by 
iGibb fc Hoarg (|200l . they find four sources (including 
IRS 1) with bipolar morphology. All appear to be domi- 
nated by wind ionized emission, although not necessarily 
jet driv en winds. K 3 -50 A, however, at a distance of 
8.7 kpc (|De Free et al.|[T994) . has a lumin osity of a mid- 
O sta r, drives an ionized bipolar outflow (jPe Free et al.l 
Il994[ ). and has a spectra l index of 0.5 in the the fre quency 
range 5 - 15 GHz (Jaff e fc Martm-Pintadol 119990 . K3- 
50 A is therefore another example of an 0-star, where the 
free-free emission appears to be coming from an ionized 
jet. Detailed studies of broad RRL objects will undoubt- 
edly discover more examples of jet-ionized H II regions. 
Such objects are likely to drive molecular outflows, excite 
masers and show evidence for strong accretion. 

To summarize: We have shown that the emission from 
IRS 1 is completely dominated by a colhmated ionized 
jet. The jet scenario also readily explains why the free- 
free emission is variable, because the accretion rate from 
the surrounding clumpy molecular cloud will vary with 
time. It also explains why IRS 1 is an extreme broad 
RRL object. Since most Hyper Compact H II regions 
are broad RRL objects, and show rising free- free emission 
with similar spectral index to IRS 1, other sources now 
classified as Hyper Compact H II regions, may also be 
similar to IRS 1. 
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